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We discuss the extension of the antenna subtraction method to include two hadrons in the initial 
state (initial-initial antennae) at next-to-next-to-leading order . We sketch the construction of the 
subtraction terms and the required phase space transformations . We discuss the integration of the 
subtraction terms in detail . 
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1. Introduction 

Jet observables are an important tool for precision studies due to their large production cross 
sections at high energy colliders. Calculating perturbative higher order corrections to jet production 
cross sections requires a systematic procedure to extract infrared singularities from real radiation 
contributions. The latter arise if one or more final state particles become soft or collinear. At 
the next-to-leading order level (NLO), several systematic and process-independent procedures are 
available . The two main methods are phase space slicing [1,2] and subtraction based methods [3,4]. 
All subtraction methods consist in introducing terms that are subtracted from the real radiation part 
at each phase space point . These subtraction terms approximate the matrix element in all singular 
limits, and are sufficiently simple to be integrated over the corresponding phase space analytically . 
After this integration, the infrared divergences of the subtraction terms become explicit and the in- 
tegrated subtraction terms can be added to the virtual corrections yielding an infrared finite result . 
Since experimental data of jet observables are reaching an accuracy of a few percent or better, ac- 
curate precision studies must rely on theoretical predictions that have the same precision . In some 
cases, this requires corrections at the next-to-next-to-leading order (NNLO) level in perturbative 
QCD. 

NNLO calculations of observables with n jets in the final state require several ingredients: the 
two-loop correction to «-parton matrix elements, the one-loop correction to (n+ l)-parton matrix 
elements, and the tree-level [n + 2)-parton matrix elements. For most massless jet observables 
of phenomenological interest, the two-loop matrix elements have been computed some time ago, 
while the other two types of matrix elements are usually known from calculations of NLO cor- 
rections to [n + 1) jet production [5] . At the NNLO level, n jet observables receive contributions 
from the one-loop (n + l)-parton matrix elements, where one of the involved partons can become 
unresolved (soft or collinear), as well as from (n + 2)-parton matrix elements where up to two 
partons can become simultaneously soft and/or collinear . In order to determine the contribution 
to NNLO jet observables from these configurations, two-parton subtraction terms have to be con- 
structed . Several NNLO subtraction methods have been proposed in the literature [6-10] . Another 
approach used for NNLO calculations of exclusive observables is sector decomposition [11] . 

In ref. [12], an NNLO subtraction method was developed for observables with partons in 
the final state only, the antenna subtraction method. It constructs the subtraction terms from the 
so-called antenna functions. The latter describe all unresolved partonic radiation between a hard 
pair of colour-ordered partons, the radiators. The antenna functions are derived systematically 
from physical matrix elements and can be integrated over their factorized phase space. At the NLO 
level, this formalism can handle massless partons in the initial or final states [16], as well as massive 
fermions in the final state [13]. For processes with initial-state hadrons, NNLO antenna subtrac- 
tion terms have to be constructed for two different cases: only one radiator parton is in the initial 
state (initial-final antenna) or both radiator partons are in the initial state (initial-initial antenna). 
Recently, in [14, 15], NNLO initial-final antenna functions were derived and integrated over their 
factorized phase space . The case with two radiators in the initial state is however still outstanding . 



2 



Initial-initial antennae 



Radja Boughezal 



In this contribution, we discuss the derivation of NNLO initial-initial antenna functions. We 
briefly describe the construction of the subtraction terms and the required phase space transforma- 
tions and discuss how the phase space integrals for initial-initial antennae can be performed using 
multi-loop techniques . 



2. Subtraction terms for initial-initial configurations 

At NNLO, there are two types of contributions to m-jet observables that require subtraction: 
the tree-level m+ 2 parton matrix elements (where one or two partons can become unresolved), and 
the one-loop m + 1 parton matrix elements (where one parton can become unresolved). In the tree- 
level double real radiation case, we can distinguish four different types of unresolved configurations 
depending on how the unresolved partons are colour connected to the emitting hard partons (see 
ref. [12] for a detailed description of the four cases). In this contribution, we focus on the case with 
two colour-connected unresolved partons (colour-connected). This is the only case where new 
ingredients are needed, namely the four-parton initial-initial antenna functions. The unintegrated 
ones can be obtained by crossing two partons to the initial state in the corresponding final-final 
antenna functions, which can be found in [12], and have then to be integrated analytically over 
the appropriate antenna phase space . The corresponding NNLO antenna subtraction term, to be 
convoluted with the appropriate parton distribution functions for the initial state partons, for a 
configuration with the two hard emitters in the initial state (partons i and / with momenta p\ and 
P2) can be written as: 

i S. colour— connected 4/ \ 1 ajt^ (1 1 \ 

da NNLO = Li d ^'«+2 (*!>•••> k m+2 \P\,Pl) 



m+2 • () '" - 

jk ~ A i, jk A Il,K~ A l,kj A iL, J ) 

jk 

X \^ m [K\,. . . ,Kl, ■ . ■ ,K m+ 2;X\Pi,X2P2)\ 2 Jm' > (Ki,. . . ,Kl, ■ ■ ■ ,K m+ 2) 



.(2.1) 



The subtraction term in eq. (2.1) is constructed such that all unresolved limits of the four-parton 
antenna function X® - k are subtracted, so that the resulting subtraction term is active only in its 
double unresolved limits, which explains the presence of the products of three-parton antennae . 
The subtraction terms for all the other unresolved configurations can be constructed using tree- 



level three-parton antenna functions. In eq. (2T) , the tree antennae X® jk , Xf - k and Xf - k depend 



on the original momenta p\, P2, kj, , whereas the rest of the antenna functions as well as the jet 
function J and the reduced matrix elements ^ m depend on the redefined momenta from the phase 

space mapping, labeled by /, / In addition to that, the reduced matrix elements depend on 

the momentum fractions x\ and X2 , which we define later . The normalization factor jV includes 
all QCD-independent factors as well as the dependence on the renormalized QCD coupling a s , 
Y,m+2 denotes the sum over all configurations with m + 2 partons, d<& m+ 2 is the phase space for an 
(m + 2)-parton final state in d = A — 2e, and finally, S m+ 2 is a symmetry factor for identical partons 
in the final state . The antenna functions can be integrated analytically, provided we have a suitable 
factorization of the phase space . The factorization is possible through an appropriate mapping of 
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the original set of momenta . These mappings interpolate between the different soft and collinear 
limits that the subtraction term regulates . They must satisfy overall momentum conservation and 
keep the mapped momenta on the mass shell . 

A complete factorisation of the phase space into a convolution of an m particle phase space 
depending on redefined momenta only, with the phase space of partons j, k, can be achieved with 
a Lorentz boost that maps the momentum q = p\+ P2~kj — kk , with q 2 > , into the momentum 
q = x\p\ +X2P2 , where X1.2 are fixed in terms of the invariants as follows [16]: 



These last two definitions guarantee the overall momentum conservation in the mapped momenta 
and the right soft and collinear behavior . The two momentum fractions satisfy the following limits 
in double unresolved configurations: 

1. j and k soft: x\ — >■ 1, X2 — >• 1, 

2. j soft and kk = Z\p\. x\ — > 1 — zi, X2 — > 1, 

3. kj =z\p\ and k k = Z2P2' *i -> 1 — zu X2 ->■ 1 — zi, 

4. kj + kk = z\p\: x\ -> 1 —zi,X2 -> 1, 

and all the limits obtained from the ones above by the exchange of p\ with P2 and of kj with kk. 
The factorized (m + 2)-partons phase space into an m-partons phase space and an antenna phase 
space is given by: 

d®m+2(kl,- ■ ■ ,k m+ 2',P\-,P2) = d®m{Kl,- ■ ■ ■ ■ ,Kk-l,Kjc + i,. . . ,K m+ 2\X\p\,X2P2) 



where [dA:] = d d k/(2n)( d ^ 8 + (k 2 ) , and J is the Jacobian factor defined by 

J =S\2 {X\{s\2-S\j-S\k)+X2{s\2-S2j-S2k)) ■ 

The next step is to integrate the antenna functions over their factorized phase space . 

3. Calculational approach for the double real radiation case 2 — > 3 

All the initial-initial antennae have the scattering kinematics p\ + p2 — > kj + kk + q, where q 
is the momentum of the outgoing particle, for example the vector boson in a vector boson plus jet 
process . Double real radiation antenna integrals are derived from squared matrix elements and can 
be represented by forward scattering diagrams as in the following figure: 





(2.2) 



x J 8(q 2 -xix 2 si 2 )S(2(x 2 p2 -x\p\).q) 
x [dkj] [dkk] dx\ dx2 , 



(2.3) 
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The two delta functions in eq.( |2.3| ) can be represented as mass-shell conditions of fake particles 
and are shown in the previous picture as a thick solid line (representing a massive particle with 
mass M = x\X2Sn) and a dashed line (representing a massless particle) . This allows us to use 
the optical theorem to transform the initial-initial antenna phase space integrals into cut two-loop 
box integrals and, therefore, use the methods developed for multi-loop calculations [17, 18] . Up 
to 8-propagator integrals with 4 cut propagators are generated in this way . The calculation of the 
integrated antennae corresponds here to the evaluation of a reduced set of master integrals . We 
found 30 of them , obtained using integration-by-part (IBP) and Lorentz identities, following the 
Laporta algorithm . We then calculate this small set of integrals using the method of differential 
equations . The simplest master integral is the two loop box with all the internal lines cut and 
defined as follows 



^^^^ =/(x 1 ,x 2 ) = Jd d qd d kjd d k k 8 d { Pl +p2-q-k j 



h) x 



8 + [k]) 5+ {k 2 k ) S + (q 2 — M 2 ) 5(2 (x 2 p 2 -x lPl ).q) . (3.1) 

As we have discussed in section [2], the phase space integrals (and therefore the master integrals) 
have to be studied in four different regions of the phase space depending on the values of x\ and 
X2, namely: 

• x\ 7^ 1, X2 7^ 1, we refer to this region as the hard one 

• x\ = 1, X2 7^ 1, and x\ 1, x% = 1, referred to as the collinear region 

• x\ = 1, x% = 1, is the soft region . 

In the hard region , the solution of the system of differential equations yields two-dimensional gen- 
eralized harmonic polylogarithms . The e expansion is needed up to transcendentality 2 . In the 
collinear regions, additional l/e coefficients may be generated and the epsilon expansion is done 
up to transcendentality 3, whereas in the soft region additional l/e 2 coefficients may appear and 
the expansion in epsilon is pushed to transcendentality 4 . We note however that the calculation 
of the masters in the soft and collinear regions, although needed with deeper expansions in e, is 
simpler than in the hard region, and only one-dimensional harmonic polylogarithms are needed in 
the collinear regions . In the soft region, a direct calculation is possible giving closed form results 
in £ in the form of gamma functions . The boundary conditions for the differential equations are 
obtained, in most of the cases, by studying the master integrals in one of the collinear limits . Oth- 
erwise the soft limit is used . 
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In a first step towards the calculation of all the integrated initial-initial antennae for the 2—^3 
tree-level double real radiation case, we have focused on all the crossings of two partons from 
the following final-final antennae: B®(q,q' ,q' ,q), E®(q,q' ,q' ,g) and H® (q,q,q ! ,q') defined in [12] , 
where the index 4 refers to four partons . There are 13 master integrals involved in their calculation, 
and the ones without irreducible scalar products are shown in Fig. [T]. 




Figure 1: Master integrals for the phase space integration of the tree-level initial-initial B%,H$ and 
£4 type antennae at NNLO . Thick solid and dashed lines refer to the conditions on the phase 
space integral implemented as auxiliary propagators . All the internal lines are massless except for 
the thick solid line . Only the integrals without numerators are shown in this picture . 



Finally , the one-loop 2—^2 antenna functions do not present any difficulty, since the needed 
one-loop box integrals are known analytically to all orders in e for the final-final configuration 
[19]. Obtaining the initial-initial contribution using these results requires crossing two legs to the 
initial state and performing the necessary analytic continuation of the involved hypergeometric 
functions . No integrals are required in this case . 

4. Outlook 

In this contribution, we have discussed the extension of the antenna subtraction formalism to 
the initial-initial configurations , including the required phase space factorisation and mappings . 
We have focused on the 2 — >■ 3 tree-level double real radiation contribution . In a first step to- 
wards the derivation of the complete set of integrated initial-initial antennae, we considered all 
the crossings of the subset of 4-parton antennae: B®(q,q' ,q' ,q), E^(q,q' ,q' ,g) and H®(q,q,q' ,q') . 
Completing the full set of NNLO antenna functions will allow the construction of subtraction terms 
needed for the evaluation of jet observables at hadron colliders . 
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